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High  Power  Laser  and  Particle  Beams  (Qiangjiguang  Yu  Lizishu) 
VoL  5  No.  3,  Aug.  1993 


SCALING  LAWS  FOR  SATURATION  INTENSITY  AND  LENGTH 
IN  AN  UNTAPERED  WIGGLER  FREE-ELECTRON-LASER  AMPLIFIER' 

Chen  Yashen,  Jiang  Youming,  and  Yang  Zhenhua 
(Institute  of  Applied  Physics  and  Computational  Mathematics, 

P.O.  Box  8009,  Beijing,  100088) 

Abstract:  A  combination  of  dimensional  analysis  and  numerical  computation  methods  was 
used  to  derive  scaling  laws  for  saturation  intensity  and  saturation  length  in  a  millimeter  wave 
free-electron-laser  amplifier  with  an  untapered  wiggler.^ 

Key  words;  saturation,  scaling  laws,  free-electron  laser  (FEL) 

1.  Introduction 

Free-electron  laser  fundamental  mode  radiation  must  pass  through  linear  and  non-linear 
stages  to  reach  saturation.  High  saturation  intensity  (the  subscript  SS  means  saturation 
value)  is  one  of  the  most  sought-after  goals.  Saturation  intensity  is  related  to  linear  stage 
growth  rate  (abbreviated  to  growth  rate  here).  Research  shows  that  when  growth  rate  is 
high,  Pss  ^  necessarily  high,  because  is  related  to  both  growth  rate  and  saturation 
length  Zss  (the  axial  distance  that  fundamental  mode  radiation  must  travel  to  reach 
saturation).  Under  conditions  approximating  one  dimension,  the  physical  parameters  that 
determine  Pss  andzgs  of  imtapered  fi-ee-electron-laser  amplifiers  are:  incident  electron  beam 
current  /,  relativity  factor  Yq,  and  relative  energy  dispersion  AYo/Voi  wiggler  magnetic  field 
strength  and  period  fundamental  mode  laser  wavelength  X  ^  and  original  fundamental 
mode  radiation  power  P^.  Of  these,  and  X  ^  satisfy  a  resonance  relationship  -  see 


This  is  the  English  title  given  by  the  Chinese  authors.  A  more  literal  translation  of  the  Chinese  title  yields 
"Scaling  laws  for  saturation  power  of  amplifier-type  free-electron  lasers."  Note  that  the  original  Chinese  title  says 
"amplifier-type  FEL"  rather  than  "FEL  amplifier." 
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Cliangcaiishu,  which  the  authors  translate  as  "untapered,"  can  also  be  translated  as  "invariable." 
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formula  (7).  Because  the  complete  expression  of  the  dependence  of  and  on  the 
above  physical  parameters  is  a  non-linear  problem,  it  is  difficult  to  use  total  analysis 
methods  on  it  In  addition,  there  are  many  physical  parameters,  so  even  if  Ayo/yq  is  given, 
there  are  also  the  four  independent  fundamental  physical  parameters  I,  Yq.  and  Yw- 
Thus,  it  is  difficult  to  use  a  pure  numerical  fitting  method  to  provide  an  expression. 

Some  newer,  better  methods  that  use  a  combination  of  dimensional  analysis  and  numerical 
simulation  are  the  United  States’  LLNL’s^  Garrison-Wong^^^  (hereafter  referred  to  as  G-W) 
and  Kumada-Sessleri^J  methods.  These  proceed  from  the  liouville  light  field  equation 
group,  and  use  the  characteristic  intensity  and  characteristic  length  of  R.  Banifacio’s 
dimensional  analysis  method  to  nondimensionalize  groups  of  equations.  The  nondimensional 
saturation  intensity  and  saturation  length  they  obtain  are  two  nondimensional  parameters 
(which  are  combinations  of  the  above  four  physical  parameters)  and  functions  of  the 
nondimensional  parameters  that  describe  the  original  conditions  (document  [1]  assumes  that 
energy  dispersion  is  zero);  afterwards,  under  conditions  where  A. 5=  l/tm  and  10/tm,  they  use 
numerical  calculation  to  approximately  determine  these  two  functions.  We  used  the  G-W 
method  to  derive  scaling  laws  for  saturation  intensity  under  millimeter  wave  conditions,  but 
our  results  were  not  as  good  as  those  attained  in  document  [1]  under  conditions  where 
and  10/^m. 

This  article  proceeds  from  a  group  of  one-dimensional  fundamental  equations  of  single 
particle  theory.  The  most  appropriate  characteristic  intensity  and  characteristic  length  were 
chosen  to  nondimensionalize  this  group  of  equations.  Results  proved  that 
nondimensionalized  single  particle  theory  equation  groups  obtained  in  this  way  included  only 
one  nondimensional  parameter,  which  is  one  fewer  than  the  G-W  method’s  nondimensional 
Liouville  light  field  equations,  and  that  there  were  also  only  two  nondimensional  parameters 
that  described  the  original  conditions:  nondimensional  original  intensity  P^q  and  A  Yq/Yo- 
this  way,  when  Ayo/Yo  given,  nondimensional  saturation  intensity  and  saturation  length 
depend  on  just  two  nondimensional  parameters.  Based  on  these  results,  it  is  easy  to  use 
numerical  calculation  to  approximately  determine  relatively  simple  scaling  laws  for  the  Bss 
and  Zgs  untapered  PEL  amplifier  with  a  certain  wavelength.  This  article  uses  this 
method  to  give  scaling  laws  for  the  P^s  and  z^s  of  millimeter  wavelengths. 
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2,  Groups  of  nondimensional  parameters  that  determine  nondimensional  saturation 
intensity  (^ss/^so)  saturation  length  (Zss/^w) 

2,1  Nondimensionalization  of  fundamental  equation  groups 

Single  particle  theory  was  used.  At  this  point,  the  one-dimensional  fundamental  equation 
group  for  untapered  wiggler  free-electron  lasers  (ignoring  the  space  charge  effect)  was  the 
electron  energy  and  phase  equation  from  document  [3]  (here,  the  subscript  j  indicates  the 
jth  electron) 

dyi  ^  _  «|g^.  .  ,  eE 

dz  yjtnc^k^  •  e- 

^  ■v  j  ®2  ^  ,  d(P 

*  •  2,'  ™Vt>.  (2) 

i'j=(k,+  k,)z-Q),t+<p+^^. 

The  amplitude  and  phase  equations  of  the  laser  electric  field  are 


sinij/.  . 

— m —  ■  \ -  /  »  l=n.ec 

m  c'k,  y.  ® 


2ne^B  /  1 


m'c’it. 


—  •  < 


In  the  above  formulas,  is  the  reduction  amplitude  of  the  laser  electric  field,  <j)  is  the  phase 
of  the  light  field,  and  are,  respectively,  the  wave  numbers  of  the  wiggler  and  the 
fundamental  mode  laser,  and  /  is  the  current  of  the  electron  beam.  For  the  circular 
polarization  and  linear  polarization  wigglers,  a  j  is,  respectively,  1  and  0.5;  a  2  is,  respectively, 
±  1  and  0.5. 

Nondimensional  electric  field 


-  L\eB, 

1 — ^  ^ 
rime- 


and  nondimensional  length 
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are  entered.  Of  these  parameters,  Lj  and  are  the  constants  of  length  dimension.  At  this 
point,  equations  (1)  through  (4)  can  be  nondimensionalized  to 


dl 


L]k  y 


^  =Lk  --^k  (1+ 


A  =  ±  / 

dF  c  ^ 


2-/; 
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m'c^k 

2/3  I?2 


me*  y»  y. 


_  1  g  ,3  L]L,k,BiI  1  . 
dz  c  me  yj  e„ 


h. 

yj 


cos  > 


(1') 

(2') 

(3') 

(4') 


2J,  Groups  of  nondimensional  parameters  that  determine  Pss/^so  ^ssAw 


From  formulas  (1’)  through  (4’),  it  is  dear  that  when  and  the  following 

resonance  relationship  is  used, 


*,(1  + 


m'c‘kl 


)  =  2k^yl 


(7) 


and  Yj  in  (1’)  through  (4’)  appears  in  the  form  of  Yj/yo,  (!’)  through  (4’)  can  therefore  be 
written  as 


d^y.hj 


ih 

n 

Ve 


=2ji[  I  y]  +  )=F,cos^,+  A 

yy  2n  yj  ‘  ^  dz 


dT 


■—  =e  <  —  > 

dz  y, 

e.y. 


±(_^  yilA^ 
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(1") 

(2") 

(3') 

(4") 

(8) 


Eqimtion  groups  (1")  through  (4")  are  equation  groups  of  nondimensional  variables  Yj/Yo> 
tjfj,  eg,  and  4>.  They  only  contain  one  nondimensional  coefficient,  e,  and  no  longer  contain 
Y  0  alone.  The  original  value  of  y  j/y  q  is  determined  by  relativistic  energy  dispersion  A  y  g/y  q; 
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the  original  value  of  ijrj  is  a  constant  that  is  unrelated  to  the  parameters  /,  and 

y  0-  The  original  value  of  <j)  is  zero,  and  the  original  value  of  ^  is  7^.  The  formula  for 
is 


e 


lO 


(9) 


In  the  above  formula,  e^o  is  the  original  value  of 

The  above  deductions  prove  that  the  nondimensional  variables  yj/y  »|r j,  and  (|)  are  merely 
functions  ofz,  e,  Ayg/y^,  and  Q)e^.  In  particular,  ej^z)  can  be  written  as 


e,(r)=/-(2.£.  4^ 


7o  v;  K 


(10) 


Because  nondimensional  saturation  length  fulfills 


, 

IJ  '--=0 


therefore,  Zss  is  just  a  function  of  <|),  ie/mc^X^BJy  and  Ay /y^.  Note  that 


/>  = 


~  B  _  , 


S=^.(^ss)  =/■(: 


^ 

■ »  IT~ 


Vo  ’  fn<^ 


e 


yl 


The  following  theorem  has  been  proved. 

Theorem  1:  For  untapered  wiggler  free-electron  laser  amplifiers, 


p  =  ^ 

^  ine*  ‘ 


ntc^  yl  ” 


70 


^SS  ~  2  (  £ 


2V  c  ,  r 

me  y 


2  ^SO  > 


^70  J 

7, 


(11) 

(12) 


In  these  formulas,  Fj  and  are  functions  of  e,  {e/mc\kjBjy  and  Ayo/yg. 
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Theorem  1  stands  up  to  arbitrary  choice  of  /,  and  Yq. 

Theorem  2:  For  conditions  where  Ayo=0  and  Fgs  has  approximately  no  relation  to  (the 
G-W  scaling  law  was  obtained  under  these  two  assumptions),  by  taking  F’i(e)oc/''^,  the  G-W 
scaling  law  is  obtained. 

Actually,  at  this  time. 


Trfe-  ( (13) 

^  tf  ^  V  y  Q 

The  above  formula  is  the  G-W  scaling  law, 

3.  Scaling  laws  for  the  saturation  intensity  and  saturation  length  of  millimeter  wavelength 
free-electron  lasers 


According  to  Theorem  1  in  part  2.2,  using  multiple  approximating  functions  within  certain 
parameter  ranges,  for  given  Ayo/y,,  and  Pgo,  Pss  can  be  written  as 


y 


(14) 


In  this  way,  the  numerical  calculation  results  of  three  /,  Yo>  and  k^  groups  with 
characteristic  meaning  (abbreviated  to  parameter  groups)  can  be  used  to  determine  the 
constants  A,  a ,  and  B  in  formula  (14),  and  then  the  computational  results  of  other  parameter 
groups  can  be  used  to  test  the  formula. 


We  used  the  above  method  on  a  millimeter  wavelength  parameter  group  range  of 
200<://Ae  cm"  2^1000,  3.5^y^nJ15,  2.75^Yw/cnK  16.5,  0.2168iAycnK2.168,  and 

13295B,„/GS(S  1626  with  an  Aurora-1  [Shuguang-1]  parameter  group  (/=500A/cm^yo=7.85, 
A„=llcm,  B,^,=4065Gs,  and  As=0,867cm)  as  its  center,  and  used  the  calculated  results  of 
document  [3]  to  derive  the  following  scaling  laws  for  saturation  intensity  and  saturation 
length  of  untapered  millimeter  wavelength  free-electron  laser  amplifiers.  When 
Fso=0,5MW/cm^  and  Ayo/yo=0.03, 
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4 

> 


/»«/ MW  =  3.672  X  10’ (£-0.203 )(  ^CB./Gs ) 

yl 

rs/cm  =  39.4  (7.,/cm  )(£-0. 203 )(  cni)’(g,/Gs) 


(15) 

(16) 


These  two  formulas  conform  well  to  the  computational  results  of  the  55  parameter  groups 
within  the  above  parameter  group  range.  For  the  great  majority,  relative  error  was  less  than 
16%,  and  for  a  small  minority,  it  was  less  than  15%.  The  range  of  change  in  e  for  these 
parameter  groups  was  0.492iei  16.7.  When  I  was  low  and  was  low,  e  would  be  low,  and 
at  this  time  was  also  low.  In  formulas  (15)  and  (16),  e- 0.203  was  a  correction  made 
whenc  was  lower,  while  under  general  conditions,  (e-0.203)«e.  Substituting  {kJ'BJyQ) 
from  formula  (16)  into  formula  (15)  yields 

i»B/MW  =  1.473x  )”«(£- 0.203  (17) 


From  the  above  formula,  it  is  clear  that  when  and  Ayo/Aq  are  given,  what  determines 
Pss  is  e  and  ZssAw»  ^  z^s- 

According  to  document  [4],  growth  rate  G(L)  is  as  follows  under  low  gain  conditions, 


Vv  4, 


(ol=  - ^ 


m 


(18) 

(19) 


while  the  highest  growth  rate  in  high  gain  situations  is 


G(L)  =  ln-^4^ 


SO 


(19) 


Note  that 


giL)  = 


icy. 


1  1  (_£_  p  KBjlV 
8  c  ^  mc»  ^ 


yl 


but  when  L=A,. 


8  c  me  yl 


(20) 
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As  a  result, 


(21) 

Therefore,  e  is  a  physical  quantity  that  describes  the  growth  rate  where  2=  In  this  way, 
according  to  the  theorems  in  part  2,  when  P^o  and  Ayo/Yo  are  given,  P^s  is  determined  by 
and  only  determined  by  growth  rate  j  and  ZgsAwJ  the  quantitative  relationship  of  millimeter 
wavelengths  is  similar  to  formula  (17).  It  is  thus  clear  that  nondimensional  saturation  length 
has  a  large  influence  on  P^. 


When  current  I  is  expressed  in  units  of  A/cm^  and  other  physical  quantities  are  expressed 
in  Gaussian  units. 


e  = 


c  me-  yl 


=  0.2  X  lO-'® 


XlBj/ 

yl 


From  formula  (15),  when  e  >  >0.203,  the  following  simplified  formula  can  yield  P^ 

^ss  =0.7344  X  10  ’  •  fi  •  /.J’5;'^y’  =  0.7344  x  I0'M„5*V 

=  0.7244UiBirr*  “  '  ^22) 

But  when  j  is  less  than  0.5,  the  error  from  simplified  formula  (22)  may  be  very  large. 

On  the  basis  of  the  preceding  theorems  and  formula  (15),  it  is  possible  to  find  the 
dependence  of  P^s  on  P^ 

4. 36677x  10’  •  />»*-’  •  (c- 0.203 )  •  (  )-^  (23) 

Vo 

We  calculated  the  above  55  parameter  groups  when  Pso= IMW/cm^  and  as  a  result,  when 
Pso=lMW/cm^  Pss  was  1  to  1.2  times  (2°'“= 1.789)  greater  than  when  Pso=0-5MW/cm^ 
More  than  half  [of  the  Pgsl  were  more  than  1.1  times  greater,  and  where  current  I  was  less 
than  300A,  [Fss]  was  less  than  1.05  times  greater.  Of  course,  when  there  are  large  variations 
in  Fso,  the  dependence  of  Fgs  0^  Aw^^w/Yo)^so  will  change,  and  the  dependence  of  P^s  on 
■^so  will  thus  change. 

Below  are  given  comparisons  of  groups  of  numerical  calculation  results  of  Fgs  Zgs  of 
parameter  groups  with  characteristic  meaning  and  comparisons  of  calculation  results  of 
formulas  (15)  and  (16),  where  Pso=0-5MW/cm^  and  Ayo/Yo=0.03.  At  the  same  time,  these 
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numerical  calculation  results  are  used  to  explain  several  physics  problems. 


(1  )7,=?.85  .  Item  ,  5.  =  4065Gs.  ;..=0.867cm9^%i!{[  / &$ 


7  ( A/  cm* ) 

1000 

625 

500 

300 

200 

/>ss(MW/cm*) 

532.8 

311.4 

257.4 

141.2 

80.3 

A, 

3.271 

3.681 

4.004 

4.5001 

4.91 

£-0.203 

9.8 

6.049 

4.796 

2.796 

1.796 

p(i) 

'  ss 

504.61 

311.49 

246.95 

143.99 

92.51 

(  J.si.  yi) 

a. 

3.326 

3.774 

3.98 

4.548 

5.036 

Key:  (1)  Current  /  varies 

Of  the  above  values,  and  ZssAw  calculated  values  from  document  [3],  and  and 
(zssAw)^^^  are  values  given  by  formulas  (15)  and  (16). 


(2)y,«7.85,  /l.=  llcm. 

1646GS  . 

=  0.2168cmV'ksit/e!c$ 

7(A/cm*) 

625 

500 

300 

7>„(MW/cm») 

155.8 

121.5 

57.9 

z„ 

A. 

7.637 

8.182 

9.82 

£-0.203 

0.822 

0.617 

0.289 

pd) 

'  ss 

164.3 

123.3 

57.8 

)... 

7.620 

8.200 

9.996 

(3)y,=  7.85,  A.=  5.5cm.  B.= 

=  8129Gs  ,  X, 

=  0.4335cm. 

/(A/cm*) 

625 

500 

300 

P^(MW/cm*) 

197.8 

150.8 

72.17 

6.28 

6.636 

7.909 

€-0.203 

1.36 

1.047 

0.547 

pd) 

*  ss 

198.1 

152.51 

79.68 

)"• 

6.37 

6.82 

8.08 

Key:  (1)  Current  I  varies 
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The  above  three  groups  of  results  were  [obtained]  under  conditions  of  different  Xg  and  A, 
when  studying  the  dependence  of  and  z^/X^  on  e.  Clearly,  when  e  >  >0.203,  is 
truly  approximately  proportional  to  j,  and  when  c  is  smaller,  P^  is  approximately 
proportional  to  €-0.203,  We  also  used  data  from  these  three  groups  to  determine  the 
degree  of  €-0.203  in  formula  (16). 

4. /  =  625A/cm*,  y,  =  7.85,  A,=  0.867cm 


A„  ( cm ) 

16.5 

14  11 

8  5.5 

(Gs ) 

2148 

2786  4065 

6655  11820 

/'ss(MW/cm>) 

339.1 

326.1  311.4 

301.9  283.3 

K 

3.242 

3.394  3.681 

4.183  4.729 

£-0.203 

8.635 

7.503  6.049 

4.483  3.143 

pin 

^ss 

342.95 

330.35  311.5 

286.53  261.65 

3.287 

3.461  3.774 

4.166  4.741 

Keys  (1)  A^  varies 

We  used  the  preceding  group  of  results  to  determine  the  degrees  and  coefficients  of 

-^w/y  0^  formulas  (15)  and  (16).  From  this  group  of  results,  one  can  see  that  when  / 

Yo,  and  A,  do  not  change  and  X^  is  reduced,  although  5, 

increases  a  great  deal,  e  still 

decreases  a  great  deal,  so  that  even  if  Zss/A„  rises,  P^  still  falls. 

(5  )  /=  625A/cm^ ,  y,=  7.85 

A„(cm) 

11 

5.5 

2.75 

-S,(Gs) 

4065 

8129 

16260 

A^(cm) 

0.867 

0.4335 

.0.2168 

/’ss(l^W/cm*) 

311.4 

197.8 

87.9 

-iaa. 

4 

3.681 

6.18 

12.2 

£-0.203 

6.049 

1.36 

0.188 

311.5 

198.2 

77.35 

3.77 

6.37 

12.38 

Keys  (1)  and  are  reduced  proportionately 
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The  group  of  results  above  shows  that  when  /  and  Yo  do  not  change  and  and  A,*  are 
reduced  proportionately,  the  results  are  similar  to  those  from  group  4:  even  though 
increases  a  great  deal,  j  decreases  very  much,  so  that  even  if  increases,  P^s  stiU 

decreases  by  a  large  margin. 


(6 )y,  =  3.92,  2.  =  2.75cm. 5. 

=  162600$  , 

,  2,= 0.867cm  ; 

I  (A/  cm’ ) 

625 

500 

300 

/>js(MW/cmO 

148.2 

110.3 

54.19 

2, 

4.546 

4.729 

5.454 

£-0.203 

2.927 

2.301 

1.299 

n(i> 

'  SS 

150.43 

118.25 

66.78 

(_^  )"’  4.527 


4.82  5.592 


Key:  (1)  /  varies 


The  preceding  is  a  miniaturized  wiggler  parameter  group  where  stays  at  0.867  cm  and 
is  increased,  while  and  Yq  are  simultaneously  decreased.  Compared  with  the  results  from 
the  first  group,  decreased  by  a  large  margin,  because  although  ZssAw  increased,  e- 0.203 
declined  a  great  deal.  Compared  with  the  fifth  parameter  group,  when  Yq  increased  from 
3.92  to  7.85,  the  Pgs  of  /=300  increased  from  54.19  to  87.9. 


(7  )y,=  11.775,  11cm  .  5.  =  29100$  .  2.=  0.2168cm 


/  ( A/  cm* ) 

625 

500 

300 

2>ss(MW/cm’) 

208.5 

136.9 

81.84 

2. 

8.09 

8.77 

10.6 

£-0.203 

0.746 

0.556 

0.253 

pill 
'  ss 

214.01 

159.56 

72.5 

)... 

8.219 

8.870 

10.89 

Key:  (1)  /varies 

The  preceding  is  a  parameter  group  where  Yo  is  increased  while  is  decreased  to  reduce 
Xj  even  more.  Compared  with  the  results  of  the  first  group,  this  happened  because  e  fell 
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by  a  large  margin.  Although  ZgsAw  inaeased,  Pss  still  declined  a  great  deal 

Finally,  we  want  to  point  out  that  from  the  computational  results  given  above,  it  is  possible 
to  find  many  examples  where  e  is  low,  whereas  2^5 A  „  is  high,  and  thus  Fgs  is  [also]  high. 

For  example,  compare  group  number  5,  where  2.75cm,  and  group  number  2,  where 
/=300A/cm^  i 

Notei  This  was  item  19075009  of  the  National  Natural  Science  Foundation  Endowment. 

The  original  of  this  document  was  received  12/24/92,  and  the  corrected  copy  was  received 
3/26/93. 
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